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ABSTRACT 


W.L: Stevens initiated research in the general area of 
random sets by considering arcs of length D (0 <D< 1) at 
random points on the circumference of a circle of length unity 
and derived the probability that the circumference is completely 
covered. D.F. Votaw Jr. considered a somewhat Meas but non- 
circular problem and derived the probability function of the 
measure of a random linear set. D.D. Beck, using the methods of 
Votaw, derived the probability function of the measure of a random 
eee set. The results of these authors are reviewed in the 
first chapter. 

The second chapter is devoted to finding the limiting 
distribution and the non-central moments of the measure of a 
random circular set for the case D = i/@tl) where’ n is the 
number of arcs. The last section of this chapter contains discussions 
of plots and tables of the probability function of the measure of 


a random circular set for various values of n and D. 
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CHAPTER I 


INTRODUCTION 


Consider ntl points Z. (i = 0,1,...,n) selected 
independently and at random from the interval (0,1); the distribu- 
tion of any Z. being the uniform distribution with distribution 


function 
CFD) F(z) =z , Ole4Z aed 


Order the ntl points in ascending order of magnitude as x, 
(i = 0,1,...,n). The totality of all cases in which two of the 
z's are equal has zero probability and such cases can be excluded 


without affecting the problem. Associate with each x, an interval 


I, where 


(x., x.+D) if. ox +Di<cd.,, 
Als al 1 
Cls2) 


tH 
i} 


+ i : 
(x, »1) u (0, x, D-1) 1£ x, +D aia ts 


where O0<D<l1 and Lien cl nis 

Let X denote the random set which is the point set sum 
of the (n+l) intervals {I} and let u(X) be its measure. The 
range of u(X) is D< u(X) < m, where m denotes the minimum of 
1 and (ntl1)D. 

D.D. Beck in his M. Sc. thesis (1968) derived‘the probability 
function of u(X) using the methods of D.F. Votaw Jr. (1946). The 


probability function, £ Oo) > of the random variable u(X) is given by 


Wiktdades” ods‘ hevhieh canes whee ae Bei at tee 


usrralow Mherenetat aoyr Sab sth a 
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6143) £&) =n. ; (n-4)? G-1)! 


n+l) mn-jt+1 
a) ets Polar 


- (1-x)I7 fx-pcjer) 23 


S 


where 


ADeceecenm( Gt) Dt gem abt 2) 4 aa MS 
M = minimum (n, [1/D]) , 


[a] = greatest integer less than or equal to a, 


and x < m. 
When m= (ntl)D, there is a finite probability that the 
Pindoneecoey Cuil condist oF nel non-overlapping intervals {I} 
and hence the distribution function, FO)» of the variable u(X) 
has a discontinuous jump or saltus point in the case when m = (n+tl1)D 
at x =m. At such a saltus point, the probability function f(s) 
of u(X), is not defined but at all other points it defines the 
probability density. The expected value of u(X) is given by 
(1.4) E(u(x)) = 1 - (1-p)™** 


and the variance is given by 


n+2 n+2 


1 
sep 2K IED) ante 2D)) ane) yeGL-D) > 
0<D <> 3 
(125) Var (u(X)) = 
2 n+2 2n+2 i 
eo (1-D)% <=(1-D) osc tie D<1. 


For D = 1/(n+1), we show in the second section of 
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Chapter II that the k non-central moment, my of the random 


variable [1-u(X)] is given by 


k Ale pau 
(1.6) om = = } Ce) fe") Ge ptekkes dg2ecivyat. . 
‘- ) 5 bi 

In the last section of Chapter II we plot and tabulate the function 
ffs) for various values of n and D. In fact, we attempted to 
find the limiting distribution of the standardized variable y(X) in 
the case D=1/(ntl1) as n _ tends to infinity, which amounts to 
increasing the number of intervals and at the same time contracting 
them. Although we were unable to determine it exactly, our results 
suggest that the limiting distribution of the variable obtained by 
standardizing u(X) is standard normal. We give the method of 
analysis we were attempting and indicate where difficulties arose, 
in the second chapter. It is unfortunate that we have been unable 
to complete the analysis of the problem, but the partial analysis 
gives strong evidence to suggest that the limiting distribution is 
normal. 

Research, which one might classify under the general head- 
ing of the measure of a random set, was initiated by a paper in 1939 
by W.L. Stevens. Beck, in his thesis (1968), gives a good account 
of the various papers published in this area since then. Consider 
the following problem: take ntl arcs of length D, O0<D<1, 
at random points on a circle of unit circumference. Let the 


location of each arc be represented by the point at its clockwise 


roltonut oda ogetyied och 1ofy oY it 


) piveaeete sv gee oh 7h ERR Om jar worm 3 


or suo. emma ah dalleltet Sa sh Alyy 2 TA 
co Sekigido ofdatyar ate te ana \ a 
, Yo Heeinad ‘ont auty em hai tbe 3 oo gaistsabasne- | 
; sees wot CNN TAES, sai oreo Stig om, gptean on anata, a 
ak dics tmod! envy ee suit enue Se foes ods erie 
abel girs arsnag at 0d ya ert Git at efquor of 


at nofssatedath gndtaett es ce soaebive grorse sewvlg 


sbaot iarondy ofa qabou meen seat oem id. seoesa 
C601 ct xsgao & et hove artas av eige ow. fans B Re,etvgn ss 903 to on 
tmsGy3e Hooy s s#viy | (AOL) alman- aiab'nk .abed . canevete IW ed) © ” 
tébtemwo°” neds Gore ote Berta nt cadet ders waeneg mitzev od? Oo 
p> 590) Ge Gignakito eouer Den niles” qealtigng gnivoltor ods 
otis 462 rT sw 36 attrib ® oe BPrtdq wolner 38 


end. Choose the end of one arc arbitrarily as the origin and 
measure distances anticlockwise around the circle. If we consider 
the point set sum of all the arcs then this representation of the 
random set is equivalent to the linear representation described 
earlier. The question is what is the probability that every point 
of the circle is included in at least one of the arcs? Or, in 
other words, what is the probability that the circle is completely 
covered? Stevens (1939) solved this problem and proved that the 


required probability is given by 


nt+1 
ah 


n+1 
2 


k nt+1 


Griz) GncaiGiep)t +07. CL-2D) Seen a) ~¢ ‘ YCI=ED) a 


where k = [1/D]. 

He also calculated the frequency distribution of the number 
of gaps. There is said to be a gap after the a are if--for ia 
distance greater than D beyond the first point of this arc, there 
is no other arc. 

In the remainder of this chapter we give the results 
obtained by D.F. Votaw Jr. (1946). Using the methods given in this 
paper, Beck (1968) derived (1.3). 

Votaw considered the following non-circular representation: 
take a random sample Z, (= 2 ene n)ed cain’ values) of one 
dimensional random variable Z having distribution function F(z). 
Arrange the values in increasing order of magnitude as Xi» 


i = 1,2,...,n and consider the set X consisting of the point set 


sum of the following intervals: 
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This representation of the random set X is not equivalent 
to the linear representation described earlier, which was equivalent 
to the circular representation. The difference arises when we con- 
sider the intervals at the end-points. Let u(X) be the measure of 
X and define a random variable y as y = u(X)-D. Then the range 
Of ey. is,,. 08S cy <u; where m, denotes the minimum of 1 and 
(n-1)D. Votaw showed that if F(z) is given by (1.1) then the 
probability function fio) of y is given by 

ee 


¥, dq J ren=1) pn-1ly) ~n-j-1 
(1.9) fy) ae ay oe (-1) ( 5 Gia ( r ) 


Mayer SyeDGay}o ae 


where 
qD < y < (qtl1)D , 
q = O,1,..-.M) : 
M, = minimum (n-2,[1/D]), and y< m, - 
He also showed that the expected value of y is given by 
_ (n-1) re eee 
(1.10) NAC) ev erapy he GED) eerie. 


Rewriting (1.9) we get 
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with the same conditions as in (1.9). Comparing (1.3) and (1.11), 
we see a lot of resemblance and the differences may be attributed 
to the following reasons: 
(i) the two representations of the random set are not equi- 

valent and hence we have two different random variables, 

(ii) Beck considered (ntl) intervals instead of the n 
considered by Votaw, and 

(iii) Beck derived the distribution of the random variable (X), 


whereas Votaw derived the distribution of u(X)-D. 


‘ 7 f F ¢ 
ae 4 
i, “.. ve 
] * _ be 
Lod 9% 
ae ee 


: iio a | 


wa oe 
pean, 9 re he eon 6. 


ee 


~tipe toa aah aie ela Sit xo note des 


nae 
ied 
: 


63. rained pobre, alee rd 
a ils io haha ee) ie | 


CHAPTER II 


MAIN RESULTS 


In this chapter, we present our attempt to find the limiting .-. 


distribution of the random variable u(X). The case D = 1/(ntl) 
is of particular interest because this amounts to reducing the size 
of each interval I and at the same time increasing the total 
number of such intervals in such a featien that it is just possible 


for the circle to be completely covered by the arcs. 


2.1 Preliminary 
When D = 1/(n+tl), the expressions (1.3), (1.4) and (1.5) 


for the probability function fox)» the expected value and the 


variance of u(X) respectively reduce to 
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Ly for all n> 1. 
Consider the limiting behaviour of the expected value and 


the variance of the random variable (XX). From (2.1.2), 


-1 


1 atl, oe ; 


0271.4) lim LL : lim {1-(1- aD 
n>o - nto 


that is, the expected value tends to a constant as n_ goes to 


infinity. However, from (2.1.3), 
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Since the limiting variance is zero, the limiting distribution of 


eae a, -1 
u(X) appears to be concentrated on the limiting mean value (l-e °). 


; 3 2 
Let us try to find the rate at which os tends to zero. 


Lemma 2.13 
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Proof: From (2.1.5), 
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2.2 Limiting Distribution 
To find the limiting distribution of a sequence of 


absolutely continuous distributions it suffices to consider the 


limit of their moment generating functions {M (t)}; where 


(212.1) M(t) = E(ett), = f enn £&) dx 
R 


R being the range of wu(X). Thus 
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Substituting the expression for £9), we get 
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When 8 = 1, we must have j r= 0 and the contribution 


to S is therefore 810 (fs hy 0 dy. When & = 2, we must have 
> 3 


j =.1 with .rx.= 0,1, or j = 2 with r= 0; the contributions. to 


S from these two cases are therefore 
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Continuing in this manner, it is readily verified that the contribu- 
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i. cae oe) Oe 
eee) Gerke) ae) ty @DIGD! 
sim 
Shue 3 


Let Ap eaten 2) 
ntl1-j-r 


3 
dy = (ne i-j—r)'dz-) Then 


i eee 


foe G3 
Pe pase ae ati ote 
M(t) ~ +l a 25 (n-j)!(5-1)! r 


0 


ae (Gg ity 


fe atjtrt(atl-j-2)2)/(mt1) (, _ jtrH(atinjor)2, 
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Caen) 


j 


nt+1 nt+1 


4-1 : es 
ety/(atl)  _ )" (ee REI aS 


so that (0 <2 <Si@and =) [@atl)-j-rlz + jrr, 


j-1 


. Pedder de Segs ce 


nt+1 


(ibys 

j=l r=0 
~ Ae 

n a e “ap? 


ntl, m—j+l) _t(jt+r)/(nt1) 
(n-j)1G-1)! ( j ee 


ie Cea dz 


Oa 1 ln Tae ae 
K | =. UO 


if “my | 


WY oa oO mt Sen OTS ae 4 


Pi s-vp 


= at 


ae ere ern ome |. " 
ee ae i: 


a3 


cs ya a ( ») (n-j) (j-1) f j (n+1-7) : i. (n+1-j-r) : 


F ' a5 

2 eh GILOHD (ptt j-r) (1 - TEE" 
je oO ig J**), ! P 
° f e oe zd Beye dz 


00 
ai) y e En.) Ger) et (itr) /(mt1) (y a Cha 
eae 
1 tr -+5z ' ; 
° | e ntl 20-4 (enon: dz 
0 
Let (l1-z) =q, dz= -dq. Then 
(2.2.4) M(t) =n iy cpt GH acd) t/t) QQ _ ity? 
oe n eae cee j-Lee OT e 4 7 ntl 
j=1 r=0 
2 ta-Tbad og ag 
| e (1-q)" 7 q?~ dq 
0 
uQO-y 
Consider the standardized variable ——~—— with moment generating 
n 
function 
((u(X)-n_)/o_)t = t/o u(X) t/o 
v(t) = Ele BT ait ese la™ete er  P 
-u_t/o 
bs fen 
=e M_(t/o_) 


From (2.2.4) 


(jon) Ut / Cea rl ey ee 
ae up) G0 a 


yo) =ne 


0 4=1 
-(tq/o_)+(tqj/(o_(mt1))) n-j oe tqr) C(ntl) 02) 
<i ~ S } Cyt Se qr n of 
a0 3 
aca 
( a aD dq 
Let 
(1 - 7 iGo ey L. a, (n) 3 so that 

(222.5) 


For y = jtr, this yields 


: = byes © 
(hs i" = @ itr) i) a, (n) Gane 
k=0 
and, since 
k 
Gtr * — (ee, A 
a lees 


we have 


: -1 ‘ = k ; * 
(1 2 itr)” r. e7 (itr) a, (a) a (e(itr)P) ; 
k=0 dp p=0 


Substituting this into the expression for v(t); we get 


x +r) n-1 
This method of dealing with the term (1 au RES was 


suggested by Dr. J.R. McGregor. 
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Jacl? oe 


I 
onl “e. ri wi hi p's I hazisom eke” 


road A 2 0d bad mpage 
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— . k 
(2.2.6) y(t) = ne a(n) ¥i(t,0) , 
where 
k 
¥<(t,0) = = ade : 
dp p=0 


and 


Ce 2ihh) ¥ ft »P) = ne 


t(1-y d/o, [ -tq/o, 
e 


0 


Ga aVee oye quan fui (athe) 
hse eo 


Wr tqr/ (ott) 6.) = Ger) + p (itr) 
Vel) alate 5 dq 
r=0 


_ ie gaan f° -tq/o, 


eo 2 Setaiaig™) qi? 
0 j=l 

Eqy Clark l)ic )etp—1 4 

e a } 


n-j ee tq/((nt1)o_) + p-l x 
Jere * Gea) He z } aq 


r=0 


t(-u d/o, i -tq/o, 


ets e 


0 


n eo ee = 
va) GiGie-o" J git 25 (a-8)®5 ag 
j=l 


where 


tq/((nt1)o_) + p-1 
(25298) Bgexe . 


Thus 


at 


( A al ant. beik 
tee 2@ wo 


ra 
Cretda +Cerky~ 6 genes per 4 ie 


tiie Bi i 


n aaa I. es 
Ft. ea) “a 7 lhe i. 


a +t rer “h 
Pe eal ts et 


ev » prt yietia 
6 


—-F 


ay 


1 


a 
«ial gay Ma ti igi ‘. 


sigh Set SONY = 


q : « ¥ 
% 
i . 
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t(1-n)/o, f -tq/o,, 


e B 
0 


(2.2.9) ¥ (tsp) = ne 


MHI) jaa = 
% ) (574) G) ca)" L(l=a). Gian J dq ° 
saz 


Consider 
a eat a 
Seay ™) ¢qe)J-* { (4-4) (4-8) FS 
fats * 4 
aoe n-1 n v -l- 
aris (Peo} (alo) aaytetesg).G-p) Hate 
v=0 
® Seal nae n-1l) mn-1+1 v 
= {(1-q) @-8) } do (emlee 2) loceien | 
Let 
@=q) (ise 2 sel et oer Cia) l= a) aap 
Then 


nn 
i} 


n-1 es 1) m-14+1) ; (s-1)/2, 
1 = {(1-q) (1-8) } o GPa Ne Nine ep wers 


5d G=a) G=8)) cS iC C5 Weer ame \ a 
Z 


(ay) a v=0 


= i 
A-ay (leh) ae 
{ (1-q) ¢ By 


CY,0) 
st1 (s) 


& 
ah 


> 


ee 


where aN is the Jacobi polynomial with parameters a=1, 


8B = 0 (see-Szegs6 (1974), equation (45332). a Now 


OL 


| 53 Ce) (Ra e, 
CAGeaty (p=) : me oP 
“Auton “on a faze 


‘a pete aa | 


, r a ei 
wis of Lei 
Ls .cwmismneg 749 Gh 


at ae 
wit 46 ia geass - wv 


iy. 


g = g8tC-g) (1-8) _  g8+(1-q) (1-8) 


~ (1-q) (1-8)-48 1-q-8 ; 
so that 
chee {oe Z 20a) 
giving 
(2.2.10) S| = (I-q-8)" + eee (au arab 


Consider the following theorem of Darboux which is given in Szegé 


(1974, page 196). 


Theorem 2.2.1: Let a and 8 be arbitrary real numbers. Then 


p6%28) -a/2 -B/2: /2 


1/2) 048 
n 


(x) ~ Gel) (x+1) ies Wa! Glas G1) 


-1/2 1/2. nF1/2 


} , 


Orn) ay ah 


(x : Be (xo 1) 


where x is outside of the closed interval [-1,1]. This formula 
holds uniformly in the exterior of an arbitrary closed curve which 
encloses the segment [-1,1], in the sense that the ratio tends 
uniformly to l. 

Lemmas. onus s = eee B lies outside the interval [-1,1]. 


tq/((at1)o,) + p-1 
Proor: Recall that 6 =e par <7 d, <p l wrandeatsco, 6 


need only be differentiated with respect to p and t in a neighbour- 
hood of p= 0, t= 0, Therefore, we may take -e < t < € and 


-6 <p <6, and choose € and 6 arbitrarily small. Thus 


S\ heey NEF ayes sie i 

| ite | 

slumsoi ead? 10; ip fay 

fate eieews baaes ene 
shaed of dar idl 


(feh+ - Xevesant. ott ebihadig avid) . 


ivq? ait) Vga" Te +e P, 
8 o@ff6 bw Lp’. & =, tanpiitacsd =i doa r 


riposigtad + htt ‘fil 5 aa oem at bo ene ® 
ol idiaiealll athe =e. ao foot " 
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ale Ss a ae Ce et xe AS ee bis 
(nt1)o Syl s (n+1)o me's Lets (n+1)o Pi ate 
n n n 
or 
-1-(eq/((nt lo) +) -1+(€q/(C(mtl1)a_)+6) 
e <B<e _ 4 
-1 


Thus e < 1/2 implies that we can choose ¢« and §& so that 


On< 2 .<. 1/24 Also "since: 0.4 q < 1, there exists dg such that 


qo = 1-8, Now 
d_ ,(1-q)(1-8)+q8, _ ___28 4 q(1-q)t 
dd [ iegae ] eigen? [l1- p+ (nto 


By restricting t to an open interval (-e€,e€), € >0 of 
arbitrarily small width, the derivative is positive over the range 
of q except when the function itself is not defined at q = 1-8. 
Now at, 
a =@)- S = 
OF Clq) G8) +98 4) ana eee tl: (1-q) (1-8) +q8 _ _} 


1-q-8 1-q-8 ; 


Therefore s lies outside the interval [-1,+1] when q « (0,1). 


Thus Theorem 2.2.1 is applicable. Note that 


oy qB+1-q-8tq8-1+qtB _ _2q8 
5 1-q-8 1-q-8 


so that 
2-1 = 498g) (1-8) 
(1-q-8) 


Then for large nf, 


i. eI | | i 
b-* “Seay 5a es 
- "4 i er 
[ 


ne, stare] \eoh-t ‘on 
f oe ot a a 


_ > ae _ 
ted ce 0 Bye 3 sat = . a 


- a iii a he 


ee 


ic 6 4% ees) tena ae e ee 
sanas ade zoe oh Bae ea nvananbieing thi oa 
E45 = Boe enttoe nce 2 sae ape ear 
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p1,0) ((1-q) (1-8)+q8) _ ¢_2q8 )-1/2 ,,2(1-q) (1-8) 7! 1/2 
Bhege a \ ieqee ) a Ge aa zee [f ise } ieee ee 
Berosye eee 
(1-q-8)* 
sf sla Eee sox nigel? apes mantayad® Perr!) 
1-q-8 lea 6 
= Rage og oe eer Seay Pere (lace ey 1/2 
Ue gee a /2 ine! edn ea? a etl acing) ee) 
v2 tql/2gh2 4 (14g) M2 (1-8)2/23 
(1-q-8) 


1/2,1/2 Vie 


375 lie + (1-q)? iT /272 


Geeta ee 


: pgllegl/ale (al 


Die iyate (Care) 


2 (yg) 2/2y2m 


n-1 pee 


} 
(iene sae 


Substituting in (2.2.10) we get 


: (qhl2gl/2 4. (1g) 4/2 (gy 2/232m 


Se ee en 
of Gas ine oe 3/4 g3/4(4- q) 1/42) 1/4 
which when substituted in (2.2.9) yields 


xe ree d/o, 
C2. 2812) y nf teP) Y 


Seen wee 


( -ta/o,, LN Peay 8 1/4 
. e 3/4 1/4 8 
0 q (1-q) 


(RAL) etmat BK Wr d-\ioe ee 
oa aN Ses 


(SVB gy? oe NS hi, 
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To evaluate this integral, we use the following generalization of 


Laplace's method (see Olver, 1974, page 333): 
Theorem 2.2.2: Let 


1 x)= r exp {-xp(t)+x”/# BCt i qt ids aa. 
0 
where k is fixed and positive, and assume that: 
(i) In the interval (0,k], p'(t) is continuous and positive, 
and the real or complex functions q(t) and r(t) are continuous. 
(11) As. t°> 30+... 


7 
p(t) = p(0) +P Eich o(t rs ; 


wa-l 
p'(t) =uP th 1 + o(t 1 ) $ 
A,-l1 
A-1 1 
q(t) = Ot -t-o(t y | 3 
v a 
rt) = Rte 7 o(t ) : where P>0O, 
Pe he 0, ST > 240, See aveceO, sand Vagn 
Then 
xp (0) 
Q Ve eA Rie af 
Hix) ooh os ) {1 v of ~ \} 
lee Lael geet pv/2 (Px) */H w/ 


w = min C4475 UyTHs v4-y) and 


F.(a,83y) is the Faxen's integral 
it 


ns r p ; ' in 
7 wi I 7 rf 


Hewes 


ay 7 
‘ae 
al 
; 


| or | 
Vv & a ‘Bes 6 Pre e! 7 a Sh “e D <a «< fh 
: 7 : = De ; - * 
; i 2} 
7 i 


Gig? * iy ame ae 2 Die (igh és 


gtBcdw ale 


ba Mery citigh W A aie wD 


aint A85 Bh ett) 
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F, (a,B3y) = | exp eee OT, 
0 
(O < Re a < 1, Re-B > 0) 


In particular 


Lerh 


: Uy re tie 
FG» 33Y) a FG: 2? y) = 20 


Ly 22 
exp G y ) 


To apply this theorem to the integral in (2.2.11), note that 


(i) jim ye = 1 e 5 
noe = 
Z 
2: -l -2 a SS 
Git)> . Lim (nt+2) o = 2e°~ -'5e so that o| =— 
ae n n n 


where Ke a aes - pe and 


tq/((ntl)o_) + p-1 


(iii) B=e southat for darge on, 6’ ~-e 


Then (2.2.11) reduces to 


ey 
a Yo/k, 1 -vn tq/k_ _(p-1)/4 
1/2 


a pee ae Ce She he at) 
nea) Geet 1)! 4 


2, oe 
¥ (tsp) as ; 


2{ (n-1) 7} 


2n 
(qe tee eieg ie o(a nee yeaa 


-1 
neiP "LD /4 1 es t/k,) (e -q) 


(O51 25) = 
of (nel) nel? yil4 0 aa 


2n 
a ae 


a Cea)e 


Gene 


where 


/2 ys 


6(q,p) = 6 = a ole +: (aye (vere) ; 


| | 
- J ¥ | | | | a 
he \\ “BNertey he @ = fe,pr0 | was 
at =“ abt tala | 


Za 


Let 


i 1 
= f exp {-n(-22n ain? E(t gy qe 4 (aq /4)7 dq . 
0 n , 


p-1y1/2 9 (P-1) /2 


Note that 6(0,p) = (l-e 4 O(L, pps = and 
e(eP +») eee es + sein Also, 
da en eet) ee ie L/ 27 weep 10/2 
ad 6(q,P) = 54 e | 5 (1-q) ier.) 
so that ‘e (etry) = 0. Let y=q- le so that when gq = 0, 
y= eee and when “q'= 1, (y=. >= apa t Thus 
(Vat /k yee") 0 ieee 
I =e { f + f } 
a _p-l 70 
| = Liz 
exp {-n(-2tn @(yteP™4,p))4+n™/? E (-y)} 
n 


Cie (eee Sea 3 


Consider I and put w= -y. Then 


(/at/k,) eee ; 


yal 1/2 
ieee | exp {-np(w)tn' r(w)}q(w)dw , 


where 


ae 


. pe Em pW (p= iy* ati 


beers B\(l-a), S < ae £6 aM : an - " 


7 “a : a | n 
pee. 
oath ok aa 


au oe Gy 


Ve son 


SM Ea Cae + 


can te 


0 = » fella Gens. oe ra a. 3 at, ‘1 = 4 
: i Me oot - ’ 


ve t me! xi 
Tih as a ieee 
| | ¥ 


or ae 


b 
Ca ae ~ sbitaiadt y 


‘Weg. 


foam: i- pe rw 4 
ee et. Se ws 
gt BY 


oe 


p(w) = 


r(w) 


Compare with Theorem 2.2.2, then 


this implies 


the value of 


also p(0) = 


" 
aut, 
EH 
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=2 an {Ctw-eP ty t/2 cP 1/2 g (Pt /2 .@-D/2, 


i Ge as oe aide 


and q(w) 


Eile 
! 
hdOo[e 


and also as_ r(w) = Ww 2 
n 
R= End ay ed) somtnocan =22) Nowstol find 


piu) — pO) , P as w-> 0+ and 


W 


Pee “note that 


QO. Thus 


Vi REY ee PAOD: = Oe an pow) = lim p'(w) 


w > O+ 


Now 
p'(w) = 


and so_ p'(0 


Also 


Let: pA = 1, 


WwW w-> 0+ w w-> 0+ 2w 


=25(itw-eP 1) M4 (q Pty 1/2 rf eae Som ee 


(d4w-eP 1/2) PH 1/2 Lely 1/2, @-D?2 
) = 0. Therefore 
pw). 7 t(1eeP“Ly-W/2 4 7 @-D)/2 


dL 
iiaha Cae’ ma 


din), Te ate 


w> Ot w 


then lim q(w) =Q which implies 
w > O+ 


gw [ed D/4 oy gPhyi/a=1 


> 


es 


SNUG), SNE Sage, SN Mey 
ot, 


Pu | 


Aa Rs! 
[OM i ny Rt Hah) © ya 
= Be i 


bali oF wot .&'= af tests Qe i =vy . 


fan <= eS 


: 2 
MM NER. | 


a CTE Pe 
ge oe ts os 
‘Nitrateg 4 Peter GR ais: 9 
j wine , ae 


a EL 
i 4 pay 
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Now applying Theorem (2.2.2), 


(Yat/k)(e*-eP*) 
e 


"Aes 2 ae 1/2 
1 2(nP) = pi/ 
Similarly, 
=p. 
Citien)(e. —e- |») 
ty = © ‘ 1/2 FG. a 
2 2(nP) < pi/ 
Thus 
Sa eet ae 
eed 
(Yate (ere) g sae Be 2 
wie Taree 20 exp GP 
2(nP) 
Substituting this into (2.2.12), we get 
5 oP D/4 (Yat/k_)(e" *-eP™*) 
ve Gtp) a 
ee eerie a 
— gqll2 9 2 (P-1)/24_ Pay 1/2 2 


ti ep-L p-l 
ET CY Nee rns Os {= e (1-e )} 
oni/? oo (P 1)/4(4_.P 11/4 Ke 
2 eP-1 


© Ph -eP*) 


(Jat/k_)(e7*-eP 7) -+ ee hie 
ae 


= exp te(p)! . 


where 


ae : re) i" Ey nao 


oe eae. f oa & 
five a. , » 

) go g ehd rire 

aU Fe 


A gaity Ts Saya + 


25 


Now 
i ¥ 6toP) = e& (P) SAD) es 
ae ©) 2 
pee ie (eee 
dp 
ae | g(p): 3 
mE) ¥ 6teP) se {(g'(p))~ + 3g" (p)e"(p)tg'" (ep) } 
P 
He @ k k-2 
SE ¥ (t.P) ~ e& P26 (g'(p))* + (2k-3)(g'(p))~ “e"(p) + ...} 
P 
Now 
p-1 2 A = ¥ 
2'(p) > Pau aes ns =; {eP Peel yea? 1), 
n k 
n 
and 
p-1 2 A a 
e"(p) = - sa oe — {eP wang: WEE Se) “ 
n k 
n 
In general, we shall have 
a Cc Ie yea (leews) Rte 
— ¥,(t,0) Te {(- BE a+ 0" )) 
dp K 
Substituting this into (2.2.6) we get 
(es Nk ye ee oe var et k 
y(t) S. & a(n) e Ge Seas. bis 
n 


which by (2.2.5) is 


es 


* 


‘teks i 


ittmans ae ee | 


38 al O.u5y ovat ods an tiveba ged, 
7 - a 


a ks MO a, Se a, 


A 
ys ih tak Maas. aoe “i 7 
q » re ah GO eta) v8 alee 
’ a) 14000 Ate rv is 7 , : 
a. ae ee 
, ae an - aa - ay a 


a / 
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(7 7k eo *(-e7) -Yate ~/k_ ete! 


n-1L 
Fh e (i + La GD) 
“4 


(7 /k Jet (1-e ) ~Yate™*/k, (n-1)an(1+(Yate™*)/ (ke, (nt1))) 
=e e e 


(07 yk e (1-e"}) ~Yate™*/k_, 


= ie e 


A (n-1) Yate” Wk (nt1)) - 5 ate /(e2(nt2)7) + 0] 
. : 


:) Syater fe 


(t” ke (1-e : l 


il 2 


5S] 
et 
@ 
b 
=S 
_ 
i 


e ”/2k-]{1+0(n"")} 


-2_ 1 


"red (t” /2k° )(2e™ Seas 
=e 


Nile 
@ 


(t" /k2){e Le 


Note that Ke +2e -5e. as n goes to infinity. Hence 


it 121 (ew 23a) Cenese)t 
lim v(t) =e 


n> 


Since for every n, v(t) is the moment generating function of 
a random variable with zero mean and unit variance the same should 


be true of lim v(t). It is clear, therefore, that at some 
n>© 


stage the approximations used have been unsatisfactory since the 
result obtained above implies that the limiting distribution is 
, = =e ae 
normal with zero mean and variance {(2e “-3e ~)/(2e “-5e “)} 
instead of unity. We have not been able to refine the approximations 


to eliminate this discrepancy. It is perhaps of some interest to 


(icgberts ONC ste yet (peat! 
es . a? 


_ fe 


fies BA contr? * wie ee 
oe ae oi noi, 


me. | (a 


‘ie "ayer de Sis\" “9 tae aa may). 


£ a t : ee * a (ee = 
' a. es se ees on! wi 
ca € oy : - 
= - a o Ube cee 3 
| OP Pustees bor sis 
— eat 2 ct 
ye 


wonal ae a == ee fn duis | 
* | 


biynite sank ‘hh aie 8 ae _ ate 
smoR If aan iotexods — e¥ a 

of goats re rere inged? ayaa ait sale 

al ieee peeiets bit 3ed%" pa of 

(e a a “at raed “aby aay | 
en ae sir sitter 02 olds ysSd You 


4 densangle sghe to sqagitad aE tt, 


zt 


note that the asymptotic evaluations of integrals by modifications 

of the Laplace method, such as we have used here, often result in 
functions which are exact except for the constants involved. Thus, 
for example, if one applies the method of Steepest Descents to deter- 
mine the F-distribution the result is exact except that the factorials 
in the normalizing constant are replaced by their Stirling's 
approximations. This may provide a clue of where one should look 


to refine the analysis to produce the exact result. 


2.3 Moments of the Variable {1-y(x)} 
For D = 1/(ntl), this section contains a result about 


the non-central moments of the random variable (1-u(X)). 


Theorem 2.3.13 The oo non-central moment of the variable 


{1-u(X)} when D=1/(n+l) is given by 


k 
1 Hh aes ae 
(2.3.4) ay = ntk ) (Ga ea (1 "ij =i) 
( k ) i=1 


pads, 2 cies = gees 


Proof: 


E(1-u(x)) * 


a 


1 k 
| (1-x) f(x) dx 
1/nt+1 \ 


n /(%+1)/(nt1) k 
| (1-x) f(s) dx 


Q=1 /2/(nt1) 


ai 6 yt 


vs 


atoiteoitibom ‘wan dennipen a . : . : 
al afvrey bse ne) ou Ae 


yeh sel? tm aa see arr” 


toul Yluede efm) erat) tt at 2 /ib nr 


Afoe ee oy sakes ea 1: soviet 


¥ 


(0) v= Eh / 


toda Yivast 6 entedeo\ tteltosa, hwy | 7 fe 


Ctx jae} 


f ip ; 
‘ t} 4 
alfdgivcay acf-lo anemog, ‘tendabartoa 


hb A 
gh I Dt De 
) ya £3) pb ME 


= ) cat ald i eas ee 


* Ina, eae gael ca a. = 


legs) 


rey 


. 
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; is /{nt1)- ; cS Gi)” be (ere) 

x ne 7 SAR hy REE SLE ‘ 

R=1 2/(nt1) j=l r=0 (n-j) © (j-1) 0 J r 
kt+j-1 ine ARS Be ae \ 


ee) ntl 


Comparing with (20222) and following the same analysis as done there, 


we get 
Bi eats A ype Gilg n+1) (n-jt+1 
(2532) ae = (n+1) dh 2 (n-j)!(j-1)! le Nel r j 
ntl k+j-1 +r n-j 
. NR 2g ager flan ba 
[ Sie Ala Se 
Let 


(n+1) P ay : 
a _ ya ktirie yas e B35 
a ed, Gl ta? (ar? dy. 


and take Zz J=i-t_. go (that “0 Si ziceley and.” gyg- (ntl=j-x) dz. 


nel=j—F 


Gee jtrt@iaj-m2 kt 
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2.4 Plots and Tables 

The results of the previous section strongly suggest that the limiting 
distribution of the standardized random variable u(X) is standard 
normal when D = 1/(ntl). In this section we examine the probability 
function, fo), of the random variable u(X) by plotting it : 
for four values of n. For each value of n, three values of 
D are considered: one less than, one equal to and one greater 


than 1/(n+l). Appendix I contains the programs used for generating 


the data, Appendix II contains the data generated and Appendix III 


vis 
anne 
; ere we 


ie yy {ee iv a a i ame 
a ; : ; “ees | 
| "eee, ae: 
; aad 4 a hari > ae 
’ a), ew * 


v a . a ¥ “ ; 
Bites Set 
: 7 ; i , . Ug 


/ . | J | | 
ya2 rimit sds reqhs Tae yt a0, ‘sien: Pa Tey oe ant - a4 
pyabes seid! AKI» soe ges vost to nolisceygekh. ge 


ot 


edidedioxy ads esitaaie ov mahtB 9, Skty ut eo = 6 Sew Loser i a 
+ ‘ghhadaly xd “Oba dhe silane: ait? is (a) 3 wotee cy 

| eo gebiny words’ - gr) to a fate unas 308 “G 20 @malev teed a fy Va 
TSIN9TRY 200 bor 62. [sipB! Sho fil aeet 406 dherbbspew0d m8 « . ww, 
attinasivg It a 4 ail peed -naie Sta «weed 


31 


contains the plots obtained by using that data. In plots 2, 3 and 
4 the scaling used for the Y axis is 1 unit = 2.0. 

Looking at the plots in Appendix III, it is clear that 
the function is symmetric only when D = 1/(nt+l) and is skewed 
in other cases. Also, even for the small values of n_ considered, 
the function tends to be a spike as n increases indicating that 


the variance decreases very rapidly. 
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GIVEN THE 


PROGRAM I 


SUBROUTINE COMB(N2,A) 
DIMENSION A(N2,N2) 

DO °SkI=475N2 

mre 1 jy a1 

NUTT) = 

CONTINUE 

NI=N2-1 

DO 16 J=2,N1 

J1i=J+1 

po 1C I=J31,N2 

A(I,J) =A (I-7,3) +A (I-14,5-1) 
CONTINUE 

RETURN 

END 
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LAST VALUE PRINTED IS (6.€,9%.9). x 
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SUBROUTINE POINTS (N2,D,A) 
DIMENSION A(N2,N2) 

REAL M1 

INTEGER Q,2 


M=INT (1/D) 

N1=N2-1 

N=N17-1 

DELS PDUTS EM) UM=N 
MI=N4*D 

BER(4? 2GT. 1) ° M1=1 

DO 15° Q=1,§ 

X=O*D 

X2=AMIN@ (M1, (Q+1) *D) 

EF. (X «ET. X2) GOTO *5 
GOTO 715 

FPNX=9 

DO 19 J=1,0 

JQ7=QO-Jt1 

Do too R=7 307 

S=A(N1,J3+1) *A(N2,5+17) *A (N2-J, 8) ¥I5 
(Ee (Raa Cake (i) 
IF( ((X-D*(J+2-7)) «EQ. (--%)) AND. 
CN 20s 0) GOTO: 1G 

S=S* (X-D* (J+R-1) ) ** (N-J) 
FNX=FNX+S “i 

CONTINUE 

WRITE (6,790) X,FNX 
FORMAT (! © F1200;. 5%, F 12.8) 
X=X+.02 

GOTO 2 

CONTINUE 

WRITE (6,19C) B1,B2 

RETURN 
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MAIN PROGRAM STARTS HERE. 


END 


EVERY TIME N IS CHANGED. RIGHT NOW 


4¢ 


DIMENSION A(21,21) 
READ(5,10) N,D1,D2,D23 
FORMAT (13,2X,3F4, 4) 


CALL 
CALL 
CALL 
CALL 
Se 
END 


COMB(N,A) 

POINTS (N,D1,A) 
POINTS (N,D2,A) 
POINTS (N,D2,A) 
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EHIS PROGRAMM PLOTS \THEe FUNCTION FNX.THE INPUT*® 
IS VALUES OF X-AND THE CORRESPONDING VALUES * 
OF THE FURCTION FNX. DATA FOR A PARTICULAR - 
VALUE OF N AND THREE VALUES OF D IS GIVEN. * 
SUE THRES SETS CF DATA ARE SEPARATED BY A * 
VALUE. OF Set, 0.0). SUBROUTINE SHOOTER IS USED* 
GO PLOT A) snooty CURVE ‘FOR SEACH SED OF DATA. * 
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CALL PLOTS 

CALL PLOT(S.¢,4.9,-3) 

CALL AXIS (0.90 ,0.9, *FNX*? ,2,6.0,90 
CALL DX Doe ig etl gp kage tg Owe pies , 


MeAND SY AXES HAVE BEEN PLOTTED, SCAGING USED is< 
Teno patent =e sy AXE SNe fae. Ce GENTH, OF 


CHO tAR nolo G UNITS EACH. 


106 FORMAT (F12.8,5X,F12.8) 
DOe30) I=1,3 
READ(5,1°9) x1,Y1 
X1=X7*%5 
V4=Y1 72.0 
CALL, SNOOTH(X1,% 17,0) 
READ(5,709) X1,Y7 

10 PSAD (5,220) X2,Y2 
DE {(K2 6EOs, Gell) cANDe (Y2neee. Oa 
* GOTO 20 
X1=X7*5 
¥i=7¥1 72.0 
Chit SMOOTH (XK) 7Y ty =) 
X1=X2 
Y1=Y¥2 
GOTO 1¢ 
26 X1=X1*5 
V=V1/2.0 
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CALL SMOOTH (X1,Y1,7-25) 
CONTINUE 

CALL PLOT(0.0,10.9,999) 
STOP 

END 
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